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Abstract 

A formula for certain exact solutions to the modified Korteweg-de 
Vries (mKdV) equation is obtained via the inverse scattering transform 
method. The kernel of the relevant Marchenko integral equation is written 
with the help of matrix exponentials as 

n{x + y;t) = Ce-<- x+y)A e SA3t B, 

where the real matrix triplet (A, B, C) consists of a constant px p matrix 
A with eigenvalues having positive real parts, a constant p x 1 matrix B, 
and a constant 1 x p matrix C for a positive integer p. Using separation 
of variables, the Marchenko integral equation is explicitly solved yielding 
exact solutions to the mKdV equation. These solutions are constructed in 
terms of the unique solution P to the Sylvester equation AP + PA = BC 
or in terms of the unique solutions Q and N to the respective Lyapunov 
equations A ] Q + QA = C f C and AN + NA^ = SB 1 " , where the f denotes 
the matrix conjugate transpose. Two interesting examples are provided. 



1 Introduction 

Consider the focusing modified Korteweg-de Vries (mKdV) equation 

1H + u xxx + 6\u\ 2 u x = 0, (1.1) 

where the subscripts denote the appropriate partial derivates, u represents a 
real scalar function and (x, t) £ M 2 . 

The modified Korteweg-de Vries (mKdV) equation arises in applications 
to the dynamics of thin elastic rods [TOJ , phonons in anharmonic lattices [52] , 
meandering ocean jets [531, traffic congestion [T3 HOI HZ1 H2] , hyperbolic surfaces 
[24], ion acoustic solitons [18], Alfven waves in collisionless plasmas [M], slag- 
metallic bath interfaces [4], and Schottky barrier transmission lines [27 . 

In this paper we present a method to construct certain exact solutions to 



(1.1 1 that are globally analytic on the entire xt-plane and decay exponentially 
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as x — > ±00 for each fixed t € R. The method used to obtain these solutions is 
based on the Inverse Scattering Transform (1ST) [J- [3], [H] and HQ. The 1ST 



matches eq.(l.l) to the first order system of ordinary differential equations 

~ = + u(x, t) 77, (1.2a) 

ax 

~ =i\r)-u(x,i)£, (1.2b) 
ax 

known as the Zakharov-Shabat system In Section [2] we develop its direct 
and the inverse scattering theory. More precisely, representing the corresponding- 
scattering data in terms of a matrix realization [5]- [7], the Marchenko integral 
equation is separable and hence, can be solved algebraically. Its solution is easily 
related to the solution of ( |1.1[ ). 

The method used in this article has several advantages: 



1. It is generalizable to the matrix version of eq.(l.l) and to other (matrix) 
nonlinear evolution equations (see, e.g., |llj). 

2. The explicit formulas found in this paper are expressed in a concise form in 
terms of a triplet (A, B, C) where A is a real square matrix of dimension 
p, C is a real row vector and B is a real column vector. When the matrix 
order of A is very large, we have an explicit formula for the solution of 



the initial value problem for (1.1). Using computer algebra, we can "un- 
zip" the solution in terms of exponential, trigonometric, and polynomial 
functions of x and t, but this unzipped expression may take several pages! 

3. Our method easily treats nonsimple bound-state poles and the time evo- 
lution of the corresponding bound-state norming constants. 

This paper is organized as follows: In Section [2] we develop the direct and the 



inverse scattering theory for system (1.2) and describe how the 1ST allows us 
to get the solution to (1.1). In Section |3| we construct the exact solutions to 



the initial value problem ( |1.1[ ) in terms of real matrices (A, B, C) solving the 
Marchenko equation. In Section|4]we write the triplet (A, B, C) in a "canonical" 
form. Finally, in Section [5] we discuss the one-soliton and a multipole solution 
as examples. 



2 1ST method for the mKdV equation 



The mKdV equation (1.1) is solvable by the inverse scattering tranform (1ST). 



This method associates ( 1.1 ) with the first order system of ODE ( 1.2 ) where the 



coefficient u is called the potential. We suppose that u is a real scalar function 
belonging to L 1 (M). In this section we recall the basic ideas behind the 1ST and 
introduce the scattering coefficients and Marchenko integral equations. A more 
detailed exposition can be found in PQ-[3], [S] and [2"T] . 

Let us introduce the J ost functions from the right tp(X,x) and ip(X,x), the 
Jost functions from the left (j)(\,x) and <^(A, x), and the Jost matrix solutions 
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*&(X,x) and <&(A,:c) from the right and from the left as those solutions of (1.2 1 
satysfying the asymptotic conditions: 



*(A,aO = ty(A,s) V(*,aO) = 
*(\,x) = (4>(\,x) 0(X,x)) = - 



e- iXJx [I 2 + 0(1)], x^+oc, 

e- iAJa; a ; (A) +o(l), x^-co, 

e- iAJa [7 2 + o(l)], 2: -►-«>, 

e~ iXJx a r {X) + o(l), x->+oo, 



where J = ^ ^ J . Here we have not written the variable t. Since this system 
is first order, we have 

$(A,x) = *(A,.x)a r (A), *(A,£c) = $(A, x)a l {X), 

where a r (A) and a; (A) are called transmission coefficient matrices. It is easy 
to prove that a r (A) and a; (A), for all A £ I, are unitary matrices with unit 
determinant, one being the inverse of the other (see, e.g., [101 02 |25]). It is 
convenient to write these matrices in the form 

lK 1 \a l3 (X) a ;4 (A)y' rV ' \a r3 (X) a r4 (X) 

where ajfe(A) and a r k(X) for k = 1,2,3,4 are real scalar functions. The Jost 
functions <f>(\,x) and V J (A, x) [TU1 GH [S3] are, for each x G R, analytic in C + , 
continuous in A G C+ and tend to 1 as |A| —¥ 00 from within C + , where we 
use C + (C~) to denote the upper (lower) complex open half planes, and = 
C ± U M. Instead, the Jost functions <f>{\,x) and ip(X,x) are, for each x G M, 
analytic in C~ , continuous in A G C~ and tend to 1 as |A| — ¥ 00 from within 
C~ . So, it is natural to consider the 2x2 matrices of functions 

F + (X,x) = (<j>(X,x) 4>(\,xj), F-(\,x) = $(\,x) 4(\,x)). 

As a result, F + (X,x) is, for each i£R, analytic in C + and continuous in <C+. 
Instead, -FL(A, x) is, for each x G M., analytic in C~ and continuous in C _ . Using 
this information we arrive at the Riemann-Hilbert problem 

F-(\,x) = F + {X,x)JS(X)J, (2.1) 

where 



i(A, "U(A) T ; (A)J 



is called the scattering matrix. 

The direct problem can be formulated as follows: Starting from the potential 
u(x), construct the scattering matrix or, equivalently, determine one reflection 
coefficient R(X) or L(X), the poles (bound-states) A.,- of the transmission coeffi- 
cient T r (X) or T;(A) (see Theorem 3.16 of [10), and the norming constants Cj S 
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corresponding to those poles. Moreover, under the technical assumption that 
aii(X), a r i(X), au(X) and a r i(X) are invertible for each A 6 1, we can relate the 
transmission coefficients to the coefficients appearing in the scattering matrix 
5(A). More precisely, we have 

T r (A) = a r i(A)-\ R(X)= -a ;4 (A)- 1 a ;3 (A) = a r3 (A)a J . 1 (A)- 1 , 
71(A) = a M (A)- 1 , L(X) = -a rl (A)- 1 a r2 (A) = a /2 (A)a i4 (A)- 1 . 

The inverse scattering problem consists of the (re)-construction of the unique 
potential u(x) if one knows the scattering data. In this paper, following [5J|5], 
we solve this problem via the Marchenko method (see [3T1[5S]) as follows: 

a. From the scattering data |i?(A), |Aj, {cjs}™1 j |, construct the 
function 

m+n n j — 1 

n(y)^R(y)+ £ X '•.-!:' ' A " ( 2 - 2 ) 

j=ro+l s— 

where = ^ J °° R(X)e lXy d\ represents the Fourier transform of R(X). 

b. Solve the Marchenko integral equation 

K(x,y)-n(x + y)^ + dz ds K (x, z) fi(z + s) Vl{s + y) t = 0, (2.3) 

J x J X 

where y > x. 

c. Construct the potential u(x) by using the following formula: 

u(x) = -2K(x,x). (2.4) 



Having presented the direct and inverse scattering problems corresponding to 
the LODE associated to the mKdV equation, we now discuss how the 1ST allows 
us to obtain the solution to the initial value problem for (1.1 1. The following 
diagram illustrates the 1ST procedure: 



given 
u{x, 0) 



mKdV 
solution 



u(x,t) 4 



direct scattering problem 
with potential u(x.O) 



inverse scattering problem 
with time evolved scattering data 



R(X), {X 3 , {^(O)}^ 1 } 



rij— l^m+n 

j—m-\-l 



time evolution 
scattering data 



°'i 



R(X,t), {XjA^mZo^Vt+i 



To solve the initial value problem for (1.1), we use the initial condition u(x,0) 
as a potential in the system (1.2). After that we develop the direct scattering 



theory as shown above and build the scattering matrix. Successively, let the ini- 
tial scattering data evolve in time. The transmission coefficient does not change 
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in time and, as a consequence, also the bound states do not. The reflection co- 
efficient R(X) evolves according to R(X, t) — e 8A t R{\). It remains to determine 
the evolution of the norming constants. Extending our previous results on the 
nonlinear Schrodinger equation ( [SJ [TTJ [TUl H] ) to the mKdV equation, we obtain 
the following time evolution of the norming constants: 



{c jnj -i{t) ... c j0 (t)) = (c jnj -i(0) ... Cj-o(O)) 
where Aj is the matrix defined by eq. (|4.7|). Finally, we solve the inverse 



scattering problem by the Marchenko method for (j 1 . 2 p with the time evolved 

> m+n 



r r — 1 1 l 

scattering data <R(X,t), <Xj, {cj S (i)}"4o f f> replacing in (2.2) R(X) 

with R(X,t) and Cj S with Cj S (t). Then the function 

u(x,t) = -2K(x,x;t), (2.5) 

is a solution to the mKdV equation. 

3 Explicit solutions for the mKdV equation 

In this section we obtain two different but equivalent formulas yielding solutions 
to the mKdV equation. Using the expressions of the evolved reflection coefficient 



and the evolved norming constants in (2.2 1, we get 

oo rn+n ry-l 

n(y;t) = - R(Xy xt e^dX+ ^ ^ c js (t)^e iX ^ , (3.1) 

which satisfies the first order PDE 

n t (y;t)+8n yvv (y;t) = 0. (3.2) 
Now, following the approach of [51 [7] , we write the kernel fl(y) introduced 



in (2.2) in the form 

n(y)=Ce- A yB, y>0, (3.3) 
where A, B 1 C are real matrices of size p x p, p x 1 and 1 x p, respectively, for 



some integer p. Eq. (3.2) suggests us to take Q(y;t) as 

Q(y; t) = Ce- Ay e 8A3t B, y>0. (3.4) 

For reasons to be clarified later, we have some further requirements on the triplet 
(A, B, C). More precisely, 

1. Our triplets (^4, B, C) is a minimal representation of the kernel £l(y), i.e, 

-f CO +oo 

p) kerCVT- 1 = p| ker B^A^y- 1 = {0}, 

r— 1 r— 1 

and we refer the reader to [8] for many details on this subject. 
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2. All of the eigenvalues of the matrix A have positive real parts. 

Following the procedure described in the preceding section, we find explicit 
solutions of the mKdV equation solving the Marchenko integral equation (2.3), 
where the kernel is given by (3.4) and the unknown function K depends on t. 
It is immediate to calculate 



fi(l/;t)t=Bt e -^v e 8(>it)»t c t j y > Q 
By using (3.4) and (3.5), eq. (2.3) becomes 

/>oo />oo 

K{x,y;t)- LBt e -At*_ / dz ds K(x,z;t)Ce- Az+8A3t e- As BBU 

J X J X 

ie _ A t s+8(A t ) s tcft = 0) y>x 



(3.5) 



-At, 



Looking for a solution of (3.6) in the form 

K{x,y-t) = H(x,t)e- A ^ +s ^ H Ci 
and introducing the matrices Q and N as 



Q = 



ds t 



-At, 



C^Ce~ As 



N 



dr 



^ — Ar 



-Atr 



after some easy calculations we obtain 

H(x,t)T(x,t) = B^e~ Atx , 
where (denoting by I p the identity matrix of order p) 



T(x, t) = Ip + e 



-A''x+a(A f ) 3 t r)p ~2Ax+SA 3 tj^ e -~A''x 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



Substituting (3.9) into (3.7) we get 



K(x,y;t) = B^e- Afx T(x 1 ty 1 e- A> y +8 ^^ t C^ = B^F-\x, t)e~ At ^C^ 
Putting 



2A t i-8(At) 3 t 

we arrive at the solution formula 

u(x,t) = -2B t F- 1 (a;,i)C t 



F{x, t) = e ZA ' x -^ A ' rt + Qe~ 2Ax+8A t N . 



(3.11) 
(3.12) 



The matrices Q, N, F(x, t) and the scalar function u(x, t) introduced by (3. 
(3.11) and (3.12), respectively, satisfy the properties stated in the following 



Theorem 3.1 Suppose that the triplet (A, B, C) is real and is a minimal rep- 
resentation of the kernel Q(y;t), and that the eigenvalues of A have positive real 
parts. Then 
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a) The matrices Q and N are real and positive self adjoint, i.e = Q and 
7Vt = N. 

b) The matrices Q and N are simultaneously invertible. 

c) The matrix F(x,t) is invertible on the entire xt-plane. Moreover, for each 
fixed t, F(x, t)^ 1 — > as x — > ±00. 

d) The real scalar function u(x,t) satisfies ( |1.1[ ) everywhere on the xt-plane. 
Moreover, u{x, t) is analytic on the entire xt-plane and decays exponen- 
tially for each fixed t as x — > ±00. 

Proof. The proof of the items a), b), c) and of the analyticity and asymp- 
totic behaviour of the solution u(x, t) is identical to the proof of items (ii) and 
(iii) of Therorem 4.2 of [7] (taking into account also Theorem 4.3 of the same 
paper). So we can refer the reader to this paper for details. However, we can 
verify directly that our solution (3.12) satisfies cq. (1.1 1. In order to do so, we 
use a slightly different notation. In particular, let us write formula (3.12) as 

u(x-t) = -2Bt e -^-r^ 1 ( a ;,i)e- Atx e 8 ( At ) 3t C t , 



(3.13) 



where T(x, t) — I p + Q(x, t)N(x) with Q(x, t) = e 



_ P -A*x „8(At) 3 t/ 



)e~ Ax e 8A * and 



N(x) 



x Ne 



-A+ a; 



, Q and N being defined in (3.8). We will see in Section 



[4] that these two matrices are the unique solutions of the so-called Lyapunov 
equations, i.e. eqs. (4.1a). We recall the following rule: If A{x) is an invertible 
matrix of functions depending on x such that its derivative with respect to x 
exists, then 

^ (A(x)- 1 ) = -A(x)- 1 ( A{x)-\ (3.14) 

Applying this differentiation rule to the function T(x, t) and deriving ( 3.13[ ) 
with respect to t , we easily get (from now on we omit the dependence of T(x, t) 
on x and f) 



lesV^r- 1 [(^) 3 - Q(x, t)A 3 N(x)] r-V^V^'c* . (3.15) 

Now, if one applies again ( |3.14[ ) to the function T and take the x-derivative, 
after some straightforward calculations and using also (4.1a), we obtain 

(r- 1 )^ = T^A^ + A^T- 1 - 2r" 1 (^ t - QAN)^ 1 . (3.16) 

As a consequence of ( |3.16[ ), we can calculate in a direct way u x obtaining 



-At a; 



T- 1 [(At) - Q(x, t)AN(x)] T- 1 e- Aix e s ^ 3t C^ . (3.17) 
Now, computing the derivative of (3.17) and taking into account ( |3.16[ ) we find 
u xx = SB^e-^zr- 1 [(A^) 2 + QA 2 N - 2(A^ - QAN)T~\A^ - QAN)] 

■ V -l e -A^ e ^ft C \ _ 
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By very similar calculations we get 



(A^ 3 - QA 3 N - 3((A t ) 2 + QA 2 N)r~ 1 (A^ - QAN) 
- Z(A* - QAN)!- 1 ^) 2 + QA 2 N) 

+ 6(A* - QAN)r~ 1 (A^ ! - QAN)^ 1 ^ - QAN)} T^e'^'S^C* . 

(3.18) 



Finally, using (3.171 and taking into account that u is a real scalar function, we 
obtain 

2\u\ 2 u x = u x u)u + uvtux = l6B^e^ A>x T^ 1 

((A*) 2 + QA 2 N)T- 1 {A^ - QAN) + (A* - QAY)r _1 ((At) 2 + QA 2 N) 

2(A* - QAN)T~ 1 (A^ - QAN)T~ 1 (A^ - QAN)]v^ 1 e' Aix e %{A ' )H C^ , 

(3.19) 



and, at this point, it is very simple to observe that u xxx + 6|u| u x 



-u t . 



Now, we will build a different explicit formula which is equivalent to the one 



expressed by (3.12). In order to obtain this result, we first observe that f2(y, t), 



as a consequence of the realness of the triplet (A, B, C), is a real function. As 
a result, 

Sl{y; i)t = n(y; t). 



By using this relation, eq. (2.3) can be written as 

(pOO poo 
Ce~ Ax - dz I dsK(x,z;t)Ce 
J X J X 

e - Ay+ 8 iA ft B = Q ^ y>x 



Az+ SA 3 t e -As BCe -Ah 



(3.20) 



By very similar computations we get the solution of (3.20) as 



K(x,y;t) = CE-\x,t)e- A ^ v - x) B. 



where 



E(x, t)=e 



2Ax-8A 3 t 



p e -2Ax+8AHp^ p = I dse -As BCe -As 



Finally, using eq. (2.5) we obtain 

v(x;t) = -2CE- 1 (x,t)B . 



(3.21) 



(3.22) 



The next theorem shows the relationships between the matrices Q, N and P 



and the solutions formula (3.12) and (3.22) 
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Theorem 3.2 Suppose that the triplet (A, B, C) is real and is a minimal rep- 
resentation of the kernel Q(y,t), and that the eigenvalues of the matrix A have 
positive real parts. Then 

i. The following relation holds: NQ = P 2 . 

ii. The matrix P is invertible on the entire xt-plane. 

Hi. E(x, t) = F(x, ty and, as a consequence, the matrix E(x, i) is invertible 
on the entire xt-plane. 

iv. For each fixed t, E(x, t)^ 1 — > as x — > ±00 . 



v. The real scalar function v(x,t) satisfy (1.1) everywhere on the xt-plane. 
Moreover, v(x,t) is analytic on the entire xt-plane and tends to zero ex- 
ponentially for each fixed t as x — » ±00. 



vi. The explicit formulas (3.12) and (3.22) yield equivalent exact solutions to 
the mKdV eq. (1.1) everywhere on the entire xt-plane. 

We skip the proof of this theorem, because it can be constructed by rearranging 
(with inessential modifications) the proofs of Theorems 5.2, 5.4 and 5.5 of [7]. 
However, we remark that: 

• By very similar calculations to those developed in the proof of Theorem 



3.1 we can directly verify that the formula (3.22) satisfies eq. (1.1) ev- 



erywhere on the xi-plane. 



• Because u(x, t) is real and scalar, the equivalence of (3.12) and (3.22) 
follows from the relation E(x, t) = F(x, ty . 

We conclude this section with the following 

Theorem 3.3 Suppose that the triplet (A, B, C) is real and is a minimal rep- 
resentation of the kernel VL{y,t), and that the eigenvalues of the matrix A have 
positive real parts. Then the solution to the mKdV equation given in the equiv- 
alent forms (3.12) and (3.22) satisfies 



[u x (x,t)Y 



d 2 log(det E(x, t)) d 2 log(det F{x, t)) 



dx 2 



dx 2 



(3.23) 



The proof of this theorem can be found in [3] . 



4 Canonical form of the triplet (A, B, C). 

In this section we show how it is possible, without loss of generality, to choose 
the triplet (A, B, C) in "canonical form." 

To obtain this representation, we begin by finding the explicit solutions of the 



mKdV given by (3.12) and (3.22) following an "algorithmic" procedure. Starting 
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with a real matrix triplet (A, B, C) which realizes a minimal representation of 
the function £l(y; t) and such that the eigenvalues of A have positive real parts, 
we consider the following equations: 



A*Q + QA = C^C, AN + NA^ = BB\ 
AP + PA = BC . 



(4.1a) 
(4.1b) 



Equations (4. la I are the so-called Lyapunov equations, instead (4.1b I is known 
as a Sylvester equation. These equations are studied in detail in [12 where it is 
proved that, under our hypotheses on the triplet (A, B, C), they are uniquely 
solvable. We have the following 



Theorem 4.1 Suppose that the triplet [A, B, C) is real and is a minimal rep- 
resentation of the kernel Q(y,t), and that the eigenvalues of the matrix A have 
positive real parts. Then, the unique solutions Q, N of the Lyapunov equations 
and the unique solution P of the Sylvester equation are such that 

1. The matrices Q, N and P are real, and Q and N are self adjoint. 



2. 



The matrices Q and N can be expressed via (3.8), instead the matrix P 
is given by the formula (3.21). Moreover, Q and N are simultaneously 



invertible, and also the matrix P is invertible. 



A proof of this theorem can be found in [7j and [12] 



Suppose that we are able to solve (4.1a) (respectively, (4.1b)). As a con- 
sequence of Theorem 4.1 their matrix solutions are those introduced in the 
preceding section by (3.8) (respectively, ( |3.21[ )). For these reasons, knowing 



the solutions Q and N (respectively, P), we can construct, in a unique way, the 



matrix F(x, t) via eq. (3.11) (respectively, E(x, t) given by (3.21)), and, finally. 



we can write down the solution of the mKdV by (3.12) (respectively, (3.22)) 



It is natural to look for a larger class including triplets such that the solutions 
of the corresponding Lyapunov or Sylvester equations have the same properties 
as in Theorem |4.1| In fact, for every triplet in this class, we can repeat the 
procedure above introduced. For this reasons we introduce the following 



Definition 4.2 We say that the triplet {A, B, C) of size p belongs to the ad- 
missible class if the following conditions are met: 

(i) The matrices A, B, and C are real. 

(ii) The triplet (A, B, C) corresponds to the minimal realization when used in 
the right-hand side of (3.3). 

(in) None of the eigenvalues of A is purely imaginary and no two eigenvalues 
of A can occur symmetrically with respect to the imaginary axis in the complex 
X-plane. 

Definition 4.3 Two triplets (A,B,C) and (A,B,C) are called equivalent if 
they lead to the same potential u(x,t). 
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Let us now consider a triplet (A, B, C) in the admissible class. What can we 
state about the solutions of the corresponding Lyapunov or Sylvester equations 



A^Q + QA^C^C, AN + NA* = BB\ 
AP + PA = BC? 

In [7J [T2] the reader will find the proof of the following 



(4.2a) 
(4.2b) 



Theorem 4.4 If the triplet (A, B, C) belongs to the admissible class, then the 
following statements hold: 



1) Equations (4.2a) and (4.2b I are uniquely solvable 



2) The matrix solutions Q and N of (4.2a) are selfadjoint. Moreover, Q and 
N are simultaneously invertible and also P is invertible. 



3) The matrices 



E(x, t) = e 



2Ax-8A-H 



Pe 



-2Ax+8A d t 



(4.3a) 
(4.3b) 



are invertible on the entire xt-plane. 
4 ) The functions 

u{x;t) = -2&'F- 1 {x,t)& ! , v(x; t) = -2CE~ l (x, t)B . (4.4) 



yield two different but equivalent explicit solutions of (1.1). Moreover, 



they are analytic on the entire xt-plane and tend to zero exponentially for 
each fixed t as x — > ±oo . 

Then for every triplet in the admissible class we can apply the algorithmic 
procedure above. 

Now, a natural question is: Starting from a triplet (A, B, C) in the admis- 
sible class, is it possible to construct an equivalent triplet (A, B, C) such that 
the matrices A, B, C are real and are a minimal representation of the func- 
tion £l(y, t), and the eigenvalues of A have positive real parts (i.e., a triplet 
(A, B, C) of type to this considered in Section [3])? The answer to this question 
is affirmative. For the sake of space we refer the reader to [7] or [B] where a com- 
plete solution to this problem is presented. In particular, eqs. (4.7) and (4.8) 
of [7] yield the triplet (A, B, C) by starting from the triplet (A, B,C), while 
eqs. (4.10), (4.11) in the same paper explain how to construct Q, N, E, F from 
Q, N, E, F. As a consequence, the following theorem allows us to understand 



which is the "canonical way" to choose the triplet (A, B, C) in (3.4) and, con 



sequently, in the explicit formulas (3.12) and (3.22). 
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Theorem 4.5 To any admissible triplet (A,B,C), one can associate a special 
admissible triplet (A,B,C), where A has the Jordan canonical form with each 
Jordan block containing a distinct eigenvalue having a positive real part, the 
column B consists of zeros and ones, and C has real entries. More specifically, 
for some appropriate positive integer m we have 



A = 


(A, 



• 

A 2 ■ 


■ M 




B = 


(BA 

B 2 




u 


• 


A m / 




\B m ) 



C = (d c 2 



(4.5) 



where in the case of a real (positive) eigenvalue ojj of Aj the corresponding blocks 
are given by 

(4.6) 



A 4 := 







:— (Cjrij 




c j2 


Cjl) , 




/^j 


-1 


••• 





o\ 









w j 


-1 • • • 










M 








UJj ■■■ 






















Bj := 

























••• 




-1 




W 







••• 





w j7 







(4.7) 



with Aj having size nj x Uj, Bj size nj x 1, Cj size 1 X rij, and the constant 
Cj nj is nonzero. In the case of complex eigenvalues, which must appear in pairs 
as aj ± ifij with aj > 0, the corresponding blocks are given by 



A, 





Cj-.= 






-h 








Aj 


-h 








Aj 











\o 














Aj 




Tji £ji) 



\ 





-h 

Aj / 



(4.8) 



B, 





VJ 



(4.9) 



where jj S and €j S for s — 1, . . . , rij are real constants with {jj n . 
denotes the 2x2 unit matrix, each column vector Bj has 2nj components, each 
Aj has size 2nj x 2nj , and the 2x2 matrix Aj is defined as 



A, 



\-f3j a 3 ) 

Proof. The real triplet (^4, B, C) can be chosen as in Section 3 of [5] 



(4.10) 
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5 Significant examples 



Example 1: Choosing the triplet (A, B, C) as 

A=(a), B=(l), C=(c) 



where a > and ^ c <G K and solving the Sylvester eq. (4.1b), we get 



By using eq. (3.22 1, we obtain 



v{x, t) 



-2c 



c z „-2ax+8a-H 
4a 2 e 



which may be called a "single-soliton solution" to (1.1). 



Example 2: In this example we consider the case in which the transmis- 
sion coefficients have a pole of order three. More precisely, let us consider the 
following triplet 



C=(l 2 1/2) 





It is not difficult to verify that the following matrix P satisfies (4.1b) 



'1/8 7/16 5/8 \ 
P = ( 1/4 3/4 13/16 
J/2 5/4 7/8 J 



In this case it is not a good idea to unzip the solution formula (3.22) in order 
to write its analytic expression because this representation take a lot of pages! 
However, using Mathematica it is very easy to plot this solution. In the next 
figure four different graphs of u(x, t) for four fixed values of t (t = 0, t = 1/4, t — 
1/2 and t = 3/4) are given. 
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